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An easily scalable toroidal geometry presents an opportunity for creating large scale persistent currents in 
Bose-Einstein condensates, for studies of the Kibble-Zurek mechanism, and for investigations of toroidally 
trapped degenerate Fermi gases. We consider in detail the process of isentropic loading of a Bose or Fermi gas 
from a harmonic trap into the scale invariant toroidal regime that exhibits a high degree of system invariance 
when increasing the radius of the toroid. The heating involved in loading a Bose gas is evaluated analytically and 
numerically, both above and below the critical temperature. Our numerical calculations treat interactions within 
the Hartree-Fock-Bogoliubov-Popov theory. Minimal change in degeneracy is observed over a wide range of 
initial temperatures, and a regime of cooling is identified. The scale invariant property is further investigated 
analytically by studying the density of states of the system, revealing the robust nature of scale invariance in 
this trap, for both bosons and fermions. We give analytical results for a Thomas-Fermi treatment. We calculate 
the heating due to loading a spin-polarized Fermi gas and compare with analytical results for high and low 
temperature regimes. The Fermi gas is subjected to irreducible heating during loading, caused by the loss of 
one degree of freedom for thermalization. 
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I. INTRODUCTION 

A range of recent experiments with ultra-cold Bose and 
Fermi gases have begun examining finite temperature and 
non-equilibrium properties in the vicinity of the phase tran- 
sition, such as critical exponents (ill, the formation of coher- 
ence |2|], and the production of defects as the system crosses 
the phase transition ytj- As system control is improved, these 
systems may provide the first quantitative experimental veri- 
fication of the scaling predicted by the Kibble-Zurek mecha- 
nism (KZM) il . 

An important and often irreducible aspect of experiments 
is the inhomogeneous potentials used to confine the gas sam- 
ples. While the most commonly realized potentials are har- 
monic there has been increasing interest in creating more gen- 
eral shapes, and here we are primarily concerned with toroidal 
traps. Sagnac interferometry is an important application for 
large toroidal traps, where the resolution is proportional to 
the area of the interferometer [51]. A variety of schemes for 
making toroidal traps using various combinations of magnetic 
and optical techniques have been proposed la]. Traps of order 
~ 3 mm diameter have been created ^ and used as a wave- 
guide for a Bose-Einstein condensate (BEC). Smaller toroids, 
in which the gas sample extends completely across the whole 
potential, have recently been produced ^^, and dynamically 
painted optical potentials iJSl] offer a new level of flexibility for 
loading non-harmonic geometries. 

Many features of toroidal geometry have been studied, in- 



cluding topological phases ]|9I], the stability of macroscopic 
persistent currents iflpll . excitation spectra llllll . atomic phase 



interference devices |12!], vortex- vortex interactions 111311 . gen- 
eration of excitations via stirring IMI . dynamics of sonic hori- 
zons il5ll . parametric arnplification of phonons lUal . rota- 
tional current generation lUTIl . the interplay of interactions and 
rotation lllSll . giant vortices Ill9ll . and vortex signatures ll20ll . 
Ideal gas theory has recently been used ll2in to study the 
rapidly rotating Bose gas in a quartically stabilized harmonic 



trap realized at ENS ]122 

The ideal gas thermodynamics of a Bose gas in a three di- 
mensional toroidal potential was first considered in ll23ll . The 
specific form of toroidal potential examined in that work was 
the harmonic-Gaussian potential that has been produced using 
combined magneto-optical traps fy\. A regime was identified 
in which the thermodynamics of the system were independent 
of the toroid radius, a feature referred to as scale invariance. 
Scale invariance should allow the efficient loading of a de- 
generate gas into a large toroidal trap, and provide a system 
uniquely suitable for quantitatively exploring the critical re- 
gion. However, scale invariance could well be highly suscep- 
tible to additional effects beyond the ideal gas treatment. In 
Ref. 12311 certain additional effects for the Bose gas were con- 
sidered: The first order interaction shift to the critical temper- 
ature was found to be scale invariant, which motivates a more 
complete treatment of the interacting system to assess the ro- 
bustness of the scale invariant predictions, and to critically 
assess the viability of this system for experimental studies. 

In this work we present a general study of the properties of 
an ultra-cold gas in the harmonic-Gaussian toroid. We gener- 
alize the ideal thermodynamics of ll23ll to a single component 
fermionic system, and develop a meanfield treatment of the 
Bose gas to account for interaction effects. The scale invari- 
ant property is further investigated analytically by studying 
the density of states of the system, revealing the robust na- 
ture and fundamental origin of scale invariance in this trap, 
for both bosons and fermions. We also give an analytical 
Thomas-Fermi treatment. The heating involved in adiabati- 
cally loading an ideal gas into a toroidal trap is evaluated an- 
alytically in the degenerate and high temperature limits for 
both the Bose and Fermi systems. We compare the analytic 
predictions for the Bose gas loading against interacting calcu- 
lations using self-consistent Hartree-Fock-Bogoliubov-Popov 
theory. Minimal heating is observed over a wide range of ini- 
tial temperatures, and a regime of cooling is identified. We 
also compare our analytical predictions for loading a spin po- 
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FIG. 1. Harmonic-Gaussian potential shown at 2; = for various 
values of Vo. For Vb < the potential is an attractive dimple trap; 
for < Vb < Vct the potential forms diflat trap (i.e. flatter at the 
origin than the pure harmonic trap); for Vb > Vo- a toroidal trap is 
realized with the minimum potential value (Vm) occurring at some 
non-zero radius (rm)- The harmonic approximation to the toroidal 
trap is shown for reference. Results for the case ctq = 24.6/im, 
Va/ks = 50nK, and Vo/ke = {-200nK, 50nK, 400nK}. 



larized Fermi gas into an approximated toroidal trap against 
a numerical calculation of the full trapping geometry. Heat- 
ing is found to be unavoidable due to the reduced degrees of 
freedom available for thermalization in the toroid. 



II. SYSTEM 

We consider a system of atoms confined in the harmonic- 
Gaussian potential 

y(x) = ^ {uoy + Lolz^) + Voexp {-ryal), (1) 

where the axial and radial trapping frequencies are w^, w^, re- 
spectively, and r is the distance from the z-axis. This trap 
can be created using a parabolic magnetic potential, com- 
bined with a detuned laser field with Gaussian spatial profile 
transverse to the propagation {z) axis, which forms an optical 
dipole potential. 

Introducing the optical potential is a relatively simple extra 
step, allowing recent experiments to use this configuration to 
great effect ylQ]- This form of harmonic Gaussian trap has 
an advantage over other realizations of toroids; changes in the 
toroid size can be easily explored by altering the laser inten- 
sity while keeping the harmonic trap and beam width param- 
eters constant, thus accessing the scale invariant regime ll23n . 

To determine the geometry of the trap it is convenient to 
define the energy 



Va 



2 2 



0' 



(2) 



which is the potential energy of an atom at r = o-q in the 
harmonic potential. In terms of the radial harmonic oscillator 
length a-r = ^Jhjmujr it takes the form 



huJr 

~2 






(3) 



so that Va expresses the area of the Gaussian in length and 
energy units of the radial trap. There are a variety of regimes 
for the harmonic-Gaussian potential, as summarized in Fig.[T] 
Here we will only consider properties of the gas either in the 
pure harmonic trap (Vq = 0) or for Vq > Va, when a true 
toroidal potential is realized. 

The minimum of the toroidal potential has a value of Vm 
and is located at z = and r = r„,, where 



K. 



as shown in Fig.[T] 



K(l + ln(l/o/K)), 
^a^,\n{Vo/V,), 



(4) 
(5) 



The radial potential can be expanded about r 



n, givmg 
harmonic trapping frequency about the minimum of the toroid 



wt 



V2i 



UJ^pT'iii 



O-Q 



= ^2\n{Vo/V^)uJr, 



(6) 



which provides the simplest approximate description of the 
radial degree of freedom. The simplest toroidal trap is then 
given by the harmonic toroid 



Vht(x) 



'■^{r - r,r,f + ujIz^ +V„, 



(7) 



characterized by only three parameters r„i, tOz, and wt, which 
are the radial location of the minimum, and the axial and radial 
trapping frequencies, respectively (see Fig.lTJ. This trap forms 
a useful approximate description in which analytical progress 
can be made. We use it extensively in what follows, and com- 
pare the results with those calculated for the full potential. 

As shown in ll23ll . quantity of fundamental interest for har- 
monic toroid traps is the geometric mean frequency of the 
toroid, w, that takes the form 



-,5 — 



W/fwfwy, 



(8) 

where hwK = fi,^/2m7rrf„ is the energy^ scale associated 
with azimuthal motion around the toroid j23i]. This geometric 
mean frequency characterizes all ideal gas thermodynamics of 
the harmonic toroid and is a central physical parameter of the 
systems studied below. 



III. BOSONS AND FERMIONS IN A HARMONIC TOROID 

Scale invariance in the harmonic toroid geometry is a rather 
specific property that should be introduced with some care. In 
this section we give a brief review of the result of Ref. |2J]: 
that the harmonic approximation to the harmonic-Gaussian 
toroidal potential exhibits scale invariance. As a point of com- 
parison we also reiterate some basic properties of the ideal 
Bose gas in a harmonic toroid. We then provide an analytical 
Thomas-Fermi treatment for the Bose and Fermi gases. 



A. Scale Invariance and density of states 



As shown in Ref. II23I1 . the concept of generalized extensive 
volume introduced by Romero-Rochin [24] can be applied to 
the harmonic toroid. The extensive volume of the harmonic 
toroid is 



V 



-5/2 



(9) 



and thus a scaling transformation of a particular degree of 
freedom that preserves V can be defined by appropriately 
changing the other variables to compensate. 

Fortunately, the scaling of lot with respect to Vb [see (|6]l] is 
precisely cancelled by the scaling of ujk in Eq. (|8]l, giving 



OJ 
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TTTTlfTn 



2'KV„/nw' 



(10) 



The absence of Vq from this expression is the basic scale in- 
variant property of the harmonic toroid: holding all param- 
eters constant while increasing Vb leads to a scaling trans- 
formation of the trap geometry which preserves the general- 
ized extensive volume. This has an intuitive explanation: as 
the toroidal perimeter increases, the radial trapping frequency 
compensates precisely by increasing so as to preserve V. 

The density of states for the harmonic toroid may be found 
asQ 



5ht(e) 



4e3/2 



3^7r(to)5' 



(11) 



where the energy, e, is measured relative to Vm- It is then 
quite simple to show that scale-invariance holds provided we 
impose the condition 



keT < K < Vo, 



(12) 



which means that the system is in a deep toroid regime de- 
scribed entirely by the density of states (fTTT i. In this regime 
the ideal gas properties (e.g. critical temperature, peak den- 
sity, entropy) are independent of Vq, and hence of the toroidal 
radius. 

However, perhaps surprisingly, it was shown in ||23I1 that 
scale-invariance in fact holds in the general harmonic Gaus- 
sian trap over a much broader range of parameters, in which 
the harmonic toroid approximation is not generally applica- 
ble. In Appendix IB1 we provide a rigorous derivation of the 
scale invariant property of the full density of states, finding 
rigorous criteria to replace the approximate set given in (fT2l) . 
[SeeEq. (EBl- 



B. Bose gas properties 

For any harmonic toroid, a simple expression for the free 
energy of the boson system can be obtained as ll23ll 



^ = iVo(y„, -m)- 



C7+/2(e 



li{p~V^) 



;3(/3to)5/2 



(13) 



where C^('Z) = SfcLi 'n^~^z^ /k^ is the polylogarithm func- 
tion. This expression yields the following estimate for the 
critical temperature for Bose-Einstein condensation in a har- 
monic toroid 



kuTr = hw 



N 



C(5/2) 



and condensate fraction 

N " 



1 



2/5 



0.89toiV2/^ 



5/2 



(14) 



(15) 



A more extensive list of thermodynamic quantities of inter- 
est is given in Ref. ll23ll . but we immediately see the conse- 
quences of changing geometry through the exponents in these 
expressions when compared with the harmonic trap, for which 

Tc oc N^/^, and Nq/N = 1 - (T/Tcf. 



C. Thomas-Fermi description of the Bose gas 

For a harmonic toroid, where the trap is completely deter- 
mined by r™, ujz and cut, and measuring /i relative to the trap 
minimum Vm, the Thomas-Fermi particle density is 



no(x) 



Uo 



1 



Rrp 



i?f 



(16) 



where this is positive, and zero elsewhere, where i?^ = 
■\/2/I/ma^, and i?f = -^/z/I/rria^ are the Thomas-Fermi 
radii, and [/q — Anfi^as/m. is the s-wave interaction param- 
eter for scattering length a^. The chemical potential can be 
expressed in terms of ([8]) as 



Mht(iVo) = feI;h/J^iVo 




(17) 



where a = y^h/mui is the geometric mean length scale 
of the toroid. The Thomas-Fermi energy is just E{Nq) = 
2/^iit(-/Vo)^o/3. Although the density changes with r^, in 
a system where uj is invariant with ?>„, the Thomas-Fermi 
chemical potential and energy are also invariants. 



D. Thomas-Fermi description of a zero-temperature 
spin-polarized Fermi gas 

We consider a spin-polarized single component Fermi gas 
at T = 0. Within the standard semi-classical description luSll 
a single component gas is described by the non-interacting 
distribution 



/(x,p) 



1 



1 



(27rft)3 e^(pV2™+v(x)-A') + i ' 



and, for example 

dxdp 



N 



{2TThy 



/(x,p) 



oPic^t^) +1' 



(18) 



(19) 



where g{e) is the density of states. For the harmonic toroid 
density of states (fTTT i. we find the Fermi energy in the har- 
monic toroid approximation 



eF,ht = ksTp 



m 



2/5 



toiV^/s ~ 1.62toiV2/5. 



(20) 

We see that the scaling with N of the toroidal system is 
slightly stronger than the N^^^ scaling of the harmonic trap 
II26I1 . and is the same as the scaling of Tc for the ideal Bose 
gas ( fT4l i. The system energy is E{T = 0) = 5eF,ht^/7. 
Integrating over momentum, the density distribution is given 
by 



n(x) 



N 



2tt RIR^Lt 



1- 



R^ 



,-i 3/2 



where this is positive, and zero elsewhere 
is the toroidal perimeter, 



R'z 



L/T — 



Rq^ 



y/2eF,ht/mu}'^ 
Fermi momentum pp 
distribution 



= ^y2€F^ht/muJ^ 
are the Thomas-Fermi radii. Defining the 



(21) 

27rr„ 

i?f = 



y^zroepYht, we find the momentum 



15 N 



'^(p) = — rr 



1-^ 

St^Pf V Pf 



(22) 



While the momentum distribution retains the characteristic 
spherical shape of the ideal Fermi gas, the power law decay 
is weaker than the (1 — p^ /pjr)^^^ behavior of the harmonic 
trap distribution. Since e^? ht is scale invariant, the peak den- 
sity in position space, and the momentum distribution are also 
invariants. 

We can also obtain a Sommerfeld expansion of the free en- 
ergy for the degenerate regime fcsT <g; fj,. Evaluating the 
grand free energy 



T ^ knT 



deg{e)\n[l + e 



-/3(e- 



giving, for the harmonic toroid density of states 



J' 






(23) 



(24) 



which gives the low temperature expansion, to 0{Pij,) 



T 



16 



105 y 7r(to)5 



357r2 



497r^ 



24(/3/^)2 384{l3fiY 



(25) 



which we use below to evaluate isentropic loading. 



IV. BOSONS AND FERMIONS IN A 
HARMONIC-GAUSSIAN TRAP 

The properties of the ideal Bose gas in the harmonic toroid 
potential have been studied extensively in ll23ll . In this section 
we obtain a convenient representation of the grand free energy 



of the ideal gas in the general harmonic-Gaussian trap that 
can be readily evaluated for both bosons and fermions. We 
also find a useful expansion of the density of states that shows 
the harmonic-Gaussian trap is generally scale invariant, and 
reveals the necessary conditions for this property. 



A. Integral representation of the grand free energy 

The various state variables can be obtained from the free en- 
ergy, which for the case of ideal gas in the harmonic-Gaussian 
potential is given by 



^0(Kn -fi)- 



T = 






bosons 



(26) 



where /3 = l/fc^T, z ~ e'^^, uj'^ ~ ^l^z, and 



fermions 



C(z,a,6) 



c« h—\ k 

a V. .. ! T{ka)-f*{ka,kb), (27) 



k=l 



k"' 



is a generalized ^-function II23I1 . with 7*(a,x) = 
x~°'^{a,x)/r{a) and 7(0,2-) = J^ e~*t°'~^dt the incom- 
plete Gamma function. When Vq -^ 0, Q{z,pVa,PVo) re- 
duces to the usual polylogarithm Q,{z). The expressions for 
the free energy are obtained by making a semiclassical treat- 
ment of the excited levels, and for the Bose case we have ex- 
plicitly indicated the condensate occupation, A^o- Consistent 
with the semiclassical approximation we take the condensate 
mode to have an energy coinciding with the potential mini- 
mum, Vm- 

As formulated in Eq. ( l27b the Q, (2, a, h) is difficult to work 
with: (i) It is not closed when differentiated with respect to 
a and h (required for obtaining the entropy and heat capac- 
ity from F etc.), rather leading to an open hierarchy of tran- 
scendental functions, (ii) The series in Eq. ( l27l i is formally 
divergent for fermions when /i > Vm ■ 

However, we have discovered a useful integral representa- 
tion of Q,{z, a, b) in terms for the usual Bose and Fermi in- 
tegrals. Using 7(a, x) = a-° L e~^'"u°-~^du, in Eq. dZTJ ). we 
obtain 



(28) 



Expression (1281 1 has the advantage of being closed under dif- 
ferentiation, and is easily evaluated numerically. For example, 
we have from S = —dJ-/dT 



Q{z,a,b) ^a / dvQ-Aze 



S_ 
ks 



V„ 



/3 (m - V^e 



dv 



K 



,/3(M-Voe-"-y, 



v)\ 



Vav 



q{^ 



Mt^-Voe-^ -Vav) 



(29) 



The representation of Eq. ( |28] ) allows all thermodynamic 
quantities for ideal Bose and Fermi gases to be calculated for 
all regimes of the harmonic-Gaussian potential (i.e. flat trap 
and dimple regimes also), and we use it to calculate the ef- 
fects of isentropic loading in SectionlVl 



B. Generalized scale invariance beyond the harmonic toroid 
approximation 

Ideal gas results can be easily determined for the more gen- 
eral harmonic-Gaussian trap once we have an expression for 
the density of states. A procedure for systematically extending 
the density of states to higher order for the harmonic-Gaussian 
trap is presented in AppendixlA] Using this approach, we find 
the density of states is scale invariant for e < Vb — Vm and is 
given by 



5(e) =5ht(e) 



1 e 

3ov: 



2520 



(30) 



to 0{e/Va-y ■ The essence of scale invariance is apparent from 
this expression: the quantity Vq, which sets the radius param- 
eter rm of the torus, does not appear, despite the distortion of 
the trap beyond the simple harmonic toroid regime. The swift 
decay of the coefficients in the series means that for most pur- 
poses the first two terms in the expansion are sufficient. 

In Appendix IBJ we use an exact expression for the density 
of states to show that when Vb > V^ the density of states is 
independent of Vq for energies < e < Vb — Vm, thus scale 
invariance formally emerges for the full harmonic-Gaussian 
trap when 



ksT « Vo - V,r. 



(31) 



We note that for a region fc^T <C Vb — Vm to exist, V^ < Vq 
is required and that the approximate condition Eq. (fT2l i is suf- 
ficient, but not necessary for (ISTT i. Our analysis also shows 
that the density of states just above Vq — Vm is almost scale 
invariant (the derivative of the density of states with respect 
to Vb goes to zero continuously at e = Vb — Vn), and so 
the temperature restriction (ISTT i is not so severe. We thus 
confirm that the scale invariant behavior of the free energy 
identified in Ref. ll23ll arises from the density of states. We 
also show in Appendix |B] that the trap density of states is 
also scale invariant, which shows that interacting extended- 
Thomas-Fermi and Hartree-Fock-Bogoliubov-Popov (HFBP) 
equations are scale invariant in the local density approxima- 
tion, see e.g. ll27ll . 

For the Thomas-Fermi chemical potential, we can see this 
explicitly by using the first few terms in gtiap{V) given in 

I which gives, for bosons, to 0([huj/Vcr]^'^/^[as/d]No) 



ti{No)=^iMiNo)<l- 
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Similarly, for the Fermi energy 



find. 



(32) 
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0([fii^/y,]i8/5Ar6/5^ 



€f = fiRht 
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iV2/5 



8 23/^ f^y'^\4f. 



441 606/5 I v^ 



(33) 



As these only depend on the invariants a and V^, general scale 
invariance is evident. 

We have shown here that scale invariance holds quite gen- 
erally, both beyond the harmonic toroid approximation, and 
for degenerate Bose and Fermi systems. Thus it should be 
possible to load degenerate quantum gases into the scale in- 
variant regime and dilate the toroid to a macroscopic scale. 
We now assess the level of heating during adiabatic loading 
from a harmonic trap (Vb = 0) into tori of different radii. 



V. ISENTROPIC LOADING 

We consider the adiabatic loading of an ultra-cold gas from 
a harmonic trap into the toroidal trap, and in what follows 
we denote thermodynamic quantities pertaining to the initial 
harmonic trap by the subscript H and those pertaining to the 
harmonic Gaussian toroid by T. Under adiabatic loading we 
equate entropy, i.e. SHiTn^N) = STiTT,N,Vo), with Vb 
the final height of the Gaussian potential. 



A. Bose gas 

1. Analytic treatment oftJie ideal Bose gas 

For the bosonic system, the cases of above and below the 
critical temperature need to be considered separately, and the 
relevant degeneracy parameter is the reduced temperature: 
tH = Th/Tch, tT = Tt/T.t, where T,t,T,h are the crit- 
ical temperatures in the initial harmonic trap, and final har- 
monic toroid trap, respectively. 

Condensed system: The regime f ^ < 1 and ty < 1 was 
considered in 12311 where it was shown that 



t^=(^S^^}^^X"tT 



V7C(3)C(7/2); 



'-H 



1.08f 



6/5 



(34) 



Uncondensed system: For the regime tn > I and It > 
1 we can make a Boltzmann approximation. In this regime 
we have that Sn/Nks -> 4 + In C(3) + 3 Ini/f and (for the 
harmonic toroid) Sr/NkB — > | — In C(5/2) + | In t^, giving 



tT 



/ C(5/2) \ 
I C(3) J 



2/5 



pl/5^6/5 



1.28t 



,6/5 
H ' 



(35) 



with the same power as in the condensed system case. While 
this analysis breaks down very close to the critical temper- 
ature (where the harmonic system may be condensed, yet 
the toroidal system uncondensed), it provides a useful upper 



bound for the jump in relative temperature that can occur in 
the vicinity of the critical point. Taking the difference between 
dMll and (ES for tn = 1, we find Air (iff = 1) = 0.19, i.e. 
a maximum change of relative temperature of less than 20%. 



2. Meanfield treatment of the interacting Base gas 

We model interactions using a HFBP treatment, which is 
a mean-field approximation, valid in the weakly-interacting 
regime. HFBP allows for interactions between the condensate 
and non-condensate and 'gives a good approximation at all 
temperatures' 12811 . 

We make a local density approximation to the HFBP, which 
is valid when the temperature is large compared to the energy 
spacing. Our approach is similar to the harmonic oscillator 
work of Giorgini et al., 1221, except that we use the general- 
ized Thomas-Fermi approximation for the condensate density 



?io(x) = -^[/i-y(x)-2[/on(x)], 



(36) 



where this is positive, and zero elsewhere, and we include the 
quantum depletion, giving 



n(x) = 



dp r /:(x, p) 



(27rft)3 i£;(x,p 

/:(x,p)-£;(x,p 

2^(x,p) 



71be[£^(x, p)] 



(37) 



where n-Q-E{E) = {e^^ — 1)~^ and £(x, p) and £'(x, p) are 
the Hartree-Fock and HFBP energy dispersions in the local 
density approximation 



£(x,p) ^^ + y(x) - /i + 2t/obio(x) + 7i(x)], (38) 
2m 



£;(x,p) 



/:2(x,p)-[C/ono(x)]' 



(39) 



We calculate the total condensate and non-condensate num- 
bers as 



-^0 = / (ixno(x), N = (ixn(x). 



(40) 



We solve Eqs. (|36]|-(|40]| self-consistently, finding /i so that 
N = Nq + N, for fixed N,m,uJr,uiz,VQ, (To and ag. Exam- 
ples of our HFBP results are shown in Fig. |2] In Fig. |2ja) we 
see the effect of increasing Vq/Vo- on the condensate fraction 
versus temperature. Notably we see a suppression of the crit- 
ical temperature Tc with increasing Vq, and the emergence of 
scale invariance, whereby the results for Vo/Va = 5 and 10 
are visually indistinguishable. In Fig. [2b) we show a typical 
example of the condensate and non-condensate density in the 
toroidal potential. 

An important effect of interactions is the depletion of the 
condensate at zero temperature, known as the quantum deple- 
tion. While the quantum depletion is generally very small in 
harmonic traps (typically < 1% OOll ). other potentials such 
as lattices can considerably enhance this IBU 13211 . By visual 
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FIG. 2. (a) Condensate fraction versus temperature for fixed 
Va/kB ~ 94nK and harmonic trap (dashed), Vo/Va = 1 (dot- 
ted), 5, 10 and harmonic toroid (all solid and indistinguishable on 
this scale). The x indicates the results studied in subplot (b). The 
critical temperatures are 20, 13 and 10 nK respectively, (b) Radial 
density at z = for the non-condensate (solid), condensate (dashed) 
and pure Thomas Fermi (dotted) [from Eq. J16b with jJ.{No) from 
J17t evaluated for the same A^o as the finite temperature system]. 
Density profiles are shown for the case indicated with x in (a), at 
T = 0.8Tc with Va/kB = 94nK and Vo/Va = 5. Calcula- 
tions are for a system of 10^ *^Rb atoms {as = 5.77 nm), with 
[u}^,LJr] = 27r[15.3, 7.8] Hz (so V„ = 200&j). 



inspection of Fig. |2|a) it is clear that the quantum depletion is 
not significantly enhanced in the toroid. We can gain some ap- 
proximate quantitative understanding of this as follows: Using 
n(x) = (TOt/onc(x))3/2/37r2ft3 for the depletion at T = 
from II33I1 . to lowest order Uanc{x) « l^htiN) — V{x), us- 
ing the chemical potential from ([TtI i. and the harmonic -toroid 
density of states from the first term of (IA9b . we get in the har- 
monic toroid limit 



Nt= 



235/8 

157r9/8 Va 



(t») 



5/4 



0.382 



V a 



N 



5/4 



(41) 



Vo/Kt 



Numerical HFBP 
Harmonic approximations 



1 5 10 harmonic toroid 



174 110 93 93 

177 



93 
94 



TABLE I. Numerical HFBP predictions of depletion at zero temper- 
ature for the cases in Fig. 13 a). The analytical estimates for a purely 
harmonic trap and a harmonic toroid are given in the bottom row. 



compared to Nt=o = 3^/^5^/^{Nas/af/y2^/^ 
0.380(iVas/a)^/^ for the harmonic trap, where a ; 



y^h/rmo 

is the geometric mean harmonic oscillator length. We find that 
these estimates agree well with the full HFBP calculations as 
shown in Table m 



3. Application of meanfield treatment to isentropic loading of an 
interacting Bose gas 

Using our meanfield theory we consider the effect of isen- 
tropic loading on the system degeneracy. To do this we ini- 
tially determine the entropy of a harmonically trapped sample 
at temperature Th, using 



S 



dxdp 

{2TThf 



{/3£;(x,p)fiBE[^(x,p)]-ln 



1 



-/JB(x,p) 



(42) 



Then, for each toroidal configuration of interest we find the 
temperature (Tt) at which the same entropy is obtained [again 
using Eq. (I42l i1. This determines the isentropic mapping 
Th -^ Tt that characterizes how the absolute temperature 
changes during loading. However, we are more interested in 
the change in degeneracy that occurs, i.e., the change in re- 
duced temperature, and for this purpose we use the critical 
temperatures Tch and Tct self-consistently determined by the 
HFBP calculation. 

Our results for the change in degeneracy are shown in Fig. [3] 
for the case of a system of 10^ ^^Rb atoms in two different 
toroidal configurations. Here we compare the change in re- 
duced temperature that occurs as a function of the initial re- 
duced temperature in the harmonic trap. 

A characteristic feature of all curves is the distinctive steep 
middle segment that occurs near t/f ^ 1. This arises as the 
gas passes through the critical point, noting that the critical 
point entropy is lower in the toroid than the harmonic trap 
(see Tab. |II|. Thus as the initial temperature tn (and hence 
entropy) increases, a point is reached (the beginning of the 
steep region) where isentropic loading will transform a par- 
tially condensed system in the harmonic trap into a normal 
system in the toroid. The steep region then terminates when 
f jy = 1 and the harmonically trapped gas also becomes nor- 
mal. 

The ideal harmonic toroid result ( l34l i predicts i^ = 0.67 as 
the temperature for which the degeneracy is unchanged upon 
loading (i.e. tr ~ ty). For tn < t*jj Eq. ( l34l i predicts that 
loading cools the system, and for tn > t*^ heating will oc- 
cur (see dotted curve in Fig. [3]l. In general we can identify 
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FIG. 3. (Color online) Effect on the reduced temperature of adi- 
abatically loading bosons into a toroid with (a) Vo- = 20/Jaj, (b) 
Va = 200fej. Other parameters are for iV = 10^ '^^Rb atoms 
with [uJ:,,uJr] = 27r[15.3, 7.8] Hz, and Vo/K = 1 (green/light- 
grey curve), 5 (red/dark-grey curve), 10 (black curve), and the har- 
monic toroid approximation (dashed curve). Also shown are the ex- 
act non-interacting results for the harmonic Gaussian trap (dashed- 
dotted curve), the harmonic toroid result Eq. l l34t (dotted curve for 
tn < 1), and the limits Eq. OSi l (dotted curve for Ih > 1). For (b), 
Vo/Va = 5 and Vo/Va = 10 are indistinguishable on this scale. 



t*jj as being where the curves in Fig.[3]intercept the horizontal 
axis (i.e. where tT^tn ~ 0). So, by inspection of our results, 
we see that interaction effects strongly suppress the amount of 
cooling that occurs and reduces the value of t*^ significantly 
(e.g. t*j^ Ri 0.4 in the harmonic -toroid limit). We note that 
for t_ff < 1 the change in relative temperature that occurs dur- 
ing loading is generally quite small. Thus the initial harmonic 
trap (in which temperature is easy to measure) is a good ther- 
mometer of the toroid in this degenerate regime and should 
provide a convenient mechanism for accurately preparing a 
loaded toroid at target degeneracy. 

For cases where the toroid is above its critical tempera- 
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VIj/V. 





1 


5 


10 


harmonic toroid 


(a) Va = 20huj 












keT^IVa 


2.2 


1.9 


1.6 


1.6 


1.7 


S{T,)/NkB 


3.6 


3.2 


3.0 


3.0 


2.9 


kBTp/V„ 


4.2 


3.6 


3.0 


3.0 


3.1 


S{TF)/NkB 


5.8 


5.4 


5.0 


5.0 


4.7 


(b) Va = 200nu) 












ksTa/V^ 


0.22 


0.13 


0.11 


0.11 


0.11 


S(Tc)/NkB 


3.6 


3.0 


2.9 


2.9 


2.9 


kBTplV. 


0.42 0.25 


0.19 


0.19 


0.19 


S(Tp)/NkB 


5.8 


4.9 


4.8 


4.8 


4.7 



TABLE II. Numerically obtained entropy and degeneracy tempera- 
tures for the systems studied in Sec.lVl 



reflecting the reduction in available degrees of freedom in the 
toroidal system. Thus, in contrast to the Bose gas, there does 
not exist a regime in which the Fermi gas undergoes cooling 
through loading into a harmonic toroid. 

High temperature system: Using the Boltzmann approx- 
imation, we find the entropy in the high temperature regime 
Sr/Nks = 7/2 + ln(150r/8) + f hifr which, with har- 
monic trap expression, Sn/NkB = 4 + In 6 + 3 In tn, gives 
the relation 



tr 



/ 256e 
V 257r 



1/5 



6/5 



1.55t 



6/5 
H ' 



(46) 



which has the same exponent as found for Bose gas loading. 



ture the effect of interactions is much less pronounced, in 
the sense that the interacting result for the harmonic toroid is 
in good agreement with the noninteracting result in the same 
limit. However, comparing the results between Fig. |3ja)-(b) 
we see that the rate of convergence to the scale invariant result 
(i.e. as Vo/Va -> oo) varies appreciably. For the parameters in 
Fig-El a) the critical temperature is quite high (see Tab.HIli and 
scale invariance emerges only when Vb is significantly larger 
than this [see Eq. (ISTTil. For the parameters in Fig.|3lb) the 
critical temperature is lower (see Tab. |II1) and the system is 
seen to approach scale invariance behavior more rapidly with 
increasing Vq. 



B. Fermi gas 

1. Analytic treatment oftJie ideal Fermi gas for the harmonic 
toroid 



Degenerate system: 
-drJ^lfi.iu- Using Eq. 
findtoO(/3/i)~'' 



The entropy is given by 5 = 
in the low temperature limit we 



'(7rAi)3 



St 

k^ ~9I3]J (fia))5 



1 + 



7tt' 



40(/3m)^ 



(43) 



Written in terms of the Fermi temperature in the toroid, TpT 
[see Eq. (I20li1. to lowest order in Tt, we find for the harmonic 
toroid 



St 
'Nk~R 



57r2 Tt 
6 TpT 



(44) 



The expression for the harmonic trap Il34ll reads Sn/Nks — 
7:^Th/Tfh, where TpH = [QNf/'^nw/kB is the Fermi 
temperature. In terms of the reduced temperatures tn = 
Th/Tfh, and It = Tt/Tft, we find that isentropic load- 
ing gives the reduced temperature 



tT 



6tH 

5 ' 



(45) 



2. Comparison with numerical treatment of the ideal Fermi gas 

Having established the high and low temperature behavior 
of a spin polarized Fermi gas under isentropic loading into a 
harmonic toroid, in this section we numerically determine the 
heating caused by loading into the general harmonic Gaussian 
toroid using the free energy ( |26] |, making use of the represen- 
tation ( l28T l. In the appropriate limits we verify the analytical 
results given above. 

For the Fermi gas it is convenient to consider the relative 
change in reduced temperature that occurs under isentropic 
loading. Results shown in Fig.|4] are for the same atom num- 
ber {N ~ 10''') and trap configurations as the bosonic results 
in Fig. [3] We note that, in Fig.|4|a) the zero temperature limit 
is markedly different from the harmonic toroid result. Includ- 
ing the next term in the scale invariant density of states, ( l30b , 

gives the following expansion of ( |45] |. to 0[N / {¥„ / hu)^) 



tT_ 
tH 



6 1 / 



3A^ 



7 Vl250(K/fiw)^ 



2/5 



(47) 



We see that, while harmonic toroid description is not partic- 
ularly good in Fig. HJa) due to the smallness of Va, the re- 
duced temperature still approaches a scale invariant limit at 
low temperature. The heating incurred when loading into a 
system with small V^r is less than that seen for loading into a 
harmonic toroid. 

The rate of convergence to the scale invariant result with 
increasing Vq/Kt is slower in Fig. Ua) than for bosons in 
Fig. Oa). Part of this effect is only apparent, due to Fig. |4] 
showing the relative change in reduced temperature. The 
actual effect is due to the fact that, while the bosonic criti- 
cal temperature (fT4l i and the Fermi temperature ( l20l i have the 
same dependence on physical parameters, the prefactor for the 
Fermi energy is almost twice as large (see Table Hill, and also 
due to the Fermi distribution: the Fermi temperature itself is 
an important energy scale for the system, and unless 



ksTpT <Vo-Vrr, 



(48) 



the system will be have occupied levels extending to regions 
of the spectrum that are not scale invariant. 
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FIG. 4. (Color online) Relative effect on the reduced temperature 
of adiabatic loading a spin polarized Fermi gas into a harmonic- 
Gaussian toroid for (a) Va = 20/k^ and (b) V„ = 200fia'. Nu- 
merical results are shown for Vo/Va = 1 (green/light-grey curve), 5 
(red/dark-grey curve), 10 (black curve) and the harmonic-toroid limit 
(dashed curve). Also shown are the Sommerfeld (horizontal dotted 
curve, from ( I45t ) and Boltzmann (dotted curve, from l l46b ) approx- 
imations for the harmonic toroid. For (b) Vq/Vo- = 5 and 10 are 
almost indistinguishable on this scale. 



VI. DISCUSSION AND CONCLUSIONS 

A. Limitations of adiabaticity 

Any defect in the trapping geometry may have a significant 
effect when interacting with, for example, a persistent cur- 
rent. The scale invariant properties we have emphasized here 
rely crucially on having a very clean trapping geometry for 
any scaling transformations of the toroid to preserve the state 
of the system. While optical traps created with continuous 
laser light can be difficult to refine to a perfect single mode 
Gaussian state, recent work on alternative procedures to pro- 
duce lattice potentials has demonstrated a way to make very 



smooth (defect free) optical potentials with pulsed femtosec- 
ond lasers 1I35I1 . 

If the system is in the scale invariant regime, conditions for 
adiabaticity when scaling the toroid may be much less strin- 
gent. The adiabaticity criterion for scaling in time iscaic would 
normally be 2tt/ujk <C iscaio, which, for large tori, would 
impose an exceedingly long timescale, as uk ~ h/mLp', 
for toroidal circumference L. For ^^Rb, the scaling time for 
L ^ 1mm is tscaie ^ lO^s. However, to the extent that the 
potential can be engineered to be highly cylindrically sym- 
metric, a radial scaling will only couple to radial excitations 
of the gas. Thus, the criteria relaxes to 27r/wT ^ ^scaic once 
the system is in the harmonic toroid regime, and it will usually 
be the case for a large toroid that uj^^ <^ ^~k ■ f'°'^ example 
(jJk/'^t ^ 10^"' for the parameters of Fig. |2b). 



B. Limitations on toroid size 

As the toroid size is increased, the largest energy scale in 
the system becomes fnuT- In order that many radial trap states 
are occupied (where the theory we have developed is applica- 
ble), the appropriate energy scale of the system must be large 
compared to fnoT'- 

(i) For bosons at zero temperature, using (fTTI i to set 

fJ-ht{N) > hujT, we get 



2^/4 

rm. < —i^y/aaaQ 



^ 



3/2 



3/2 
(49) 



(ii) For fermions at zero temperature, using (l20t to set 
EFjit > hhJT, we get 



.™,<4('i^^V''^ao«1.14.o-iV^/^ (50) 



V2V 






(iii) In the high temperature limit we set fc^T ^ TujJt to get 

(51) 



Trr,. < -^- ~ O.lVllaQ- . 



y/2 hWr 

For the parameters of Figs. |3] and |4l the upper bounds for r^ 
are ~ 0.1 mm at T = and ~ 1 mm at Tc or Tp. By optimiz- 
ing the theory we have developed here, much larger toroids 
would be possible. 



C. Investigating KZM using scale-invariance 

The following are desirable features of the scale invariant 
toroidal trap for studies of the Kibble-Zurek effect in an ultra- 
cold Bose gas: 

a. System control: The trap should reliably produce sys- 
tems near the critical temperature and have good control over 
quenching the system through the critical point. Cooling in 
a harmonic trap and then loading using scale-invariance pro- 
vides this feature, as a target degeneracy can be produced. 
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b. Suppression of post-quench dynamics: An important 
practical issue in KZM experimental studies concerns the 
post-quench recombination dynamics of vortices: Evapora- 
tive cooling across the BEC transition in a harmonic trap is 
followed by rapid decay of the vortex density Ii36i1 , associ- 
ated with recombination of vortex anti- vortex pairs 03711 and 
the decay of lone vortices IBSll . This can greatly complicate 
measurements of defect density scaling with quench rate. The 
same process in a toroidal trap will aggregate defects into a 
macroscopic persistent current ll39ll . It may be easier in prac- 
tice to measure these energetically stabilized vortices than to 
characterize the decay of vortices in a harmonic trap. 

c. Effective one dimensional quench: A thermal gas 
quenched in a large toroid will not be quenched in the trans- 
verse degrees of freedom, because they respond much more 
rapidly than the azimuthal modes. Above Tc the relaxation 
time for mode with energy e is well estimated as ImXl 



ksT 



(52) 



where eq = corresponds to the ground state, and Hojk ^ 
hojT, fiojz will typically hold. Here 7 is the damping rate aris- 
ing from thermal cloud interactions. The system can be out of 
equilibrium azimuthally, but still in equilibrium with respect 
to transverse degrees of freedom with much higher energy. 

d. Correlation length: An essential quantity of interest 
in KZM is the correlation length ^, over which the first order 
correlation function decays. For an ideal homogeneous gas in 
thermal equilibrium at Tc < T, a semi-classical calculation 
gives im 



e 



2m\fj,\ 



(53) 



In a local density description of the trapped gas, we have 
^(x)2 = ti^/2m\fi{x.)\, with ^(x) ^ fi - V{x.) the local 
chemical potential. A scale invariant ideal gas will have a 
scale invariant chemical potential at the trap center, and hence 
the equilibrium correlation length will also be scale invariant 
there. Thus, defect formation for a range of tori with different 
perimeters could be explored for a given initial temperature 
and chemical potential of the gas (and therefore identical cor- 
relation length). For a larger toroid a larger number of correla- 
tion lengths could be frozen in during the quench (see Fig.|5]l. 
This could provide a means to make persistent currents with 
high winding number 



D. Conclusions 

In this work we have estabUshed the scale invariance prop- 
erty of a harmonic Gaussian trap, extending the largely ideal 
Bose gas treatment of ll23ll to systems of interacting bosons 
and spin polarized fermions by analyzing the density of states. 
We have assessed the level of heating incurred through adia- 
batic loading of bosons with interactions, and of spin polar- 
ized fermions. 

Specifically, we have provided a formal proof of scale 
invariance of the general harmonic-Gaussian toroidal trap. 




FIG. 5. Equipotential surfaces in for the case of a (a) haiTnonic trap, 
(b) small toroidal trap, and (c) large toroidal trap. The correlation 
length is indicated (see text). 



demonstrating that large scaling transformations can be ap- 
plied to the toroid without altering thermal equilibrium prop- 
erties due to the invariant nature of trap geometry. Focus- 
ing on adiabatic loading from a harmonic trap, we have con- 
sidered a range of possible tori determined by the Gaussian 
height Vq and the energy Va = ma;,^(To/2 which are the two 
energies defining the general trap. Analytical and numerical 
results for a single component Fermi gas, and for an interact- 
ing Bose gas are compared, allowing the relative temperature 
change under isentropic loading to be determined. 

Bose gas: The Bose gas always exhibits a regime where de- 
generacy increases slightly, occurring for initial reduced tem- 
peratures lying in the range < t// < 0.4. We also ob- 
serve that for i_y < 1 the degeneracy change that occurs using 
loading is quite small, which should allow loading of a well- 
characterized system in a harmonic trap into a toroid with the 
desired temperature. 

Fermi gas: Loading into a toroid always causes heating. 
For the special case of loading into a harmonic toroid, the 
smallest heating that may occur is 20%, reflecting the fact that 
the angular position coordinate is not available for thermaliza- 
tion in the toroidal system. Also, the Fermi energy itself sets 
an important energy scale for the onset of scale invariance, 
and that the Fermi gas shows a slower approach to the scale 
invariant regime than an equivalent Bose gas, with respect to 
increasing V^- 

We emphasize that, for finite N, the system is always scale- 
invariant as Vq/Vct — > 00, since the energy scale determining 
invariance, Vq — Vm, tends to 00. This property suggests that 
the harmonic Gaussian trap may provide a means to generate 
degenerate systems characterized by a macroscopic toroidal 
perimeter, or to make persistent currents of high winding num- 
ber. 
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Appendix A: Density of states series expansion 

Here we obtain a formal series expansion of the density of 
states in the toroidal potential, in the scale-invariant regime. 
We make use of a decomposition of the general density of 
states problem into a sequence of convolutions of densities of 
states corresponding to each additive term in the single parti- 
cle energy 113211 . We can write, for example 



g{e)=j;gK{K)gtr.p{e^K)dK, 



(Al) 



and work with the trap density of states 5trap(V^) separately 
from the kinetic density of states gK{K). We can further de- 
compose the trap density of states into radial and axial parts 

gtr.piV) = 27: f gr{Vr)g,{V ~ Vr)dVr. (A2) 
Jo 

For Vo > Va, Kn = Va[l + ln(Vb/V;)] and the effective trap 



Our main task is to find a useful expression for griV) 



(A3) 



(A4) 



1 _ ^e-«V-e 1 - 



-R''_/al 



where R± are defined to satisfy Vcfi{R±) = V, which we im- 
plicitly differentiate using (IA3l l. This requires that < i?_ < 
i?-(_, i.e. V < Vq ~ Vm, which is a necessary condition for the 
following. We proceed by expanding the potential about the 
minimum 



i=2 



r V 



(A6) 



where, for example, C2 = ■muj^/2 and C3 = ujT{3ujf — 
iij^)y/m^/Va/6. Inverting the series ll42ll . we find, to 
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FIG. 6. Density of states for the parameters of Fig. (S^b) and with 
Vo = 5Vcr'- harmonic trap (dotted), series expansion with one 
(dashed), and two (solid) terms; with two terms, the series expan- 
sion and the exact numerical result (also solid) are indistinguishable 
on this scale. 



Using Eq. (|A5]l we find, to OiV/Va 



^5/2 



gAV) 



V2 



1 V 



1 



216 V Vr 



V 



3/2' 



and ( IA2I ) gives trap density of states, to 0{V/Vcr) 
47r2^/K 



Strap (1^) 



(muj^)^/^ 



1 ^ 
IH 

12 K 



- - 

576 I 14 



(A8) 



, (A9) 



(A3) using gziVz) = ■\/2/ma;^I4- Finally, we make use of 
y/K/2/Tr'^h^ and dATl i to find the total den- 



,3/2 



gxiK) 

sity of states expansion ( l30l l. The expansion for the density of 
states is compared with exact numerical results in Fig.|6]. 



Appendix B: Proof of scale invariance 

In AppendixlAl we derived the first few terms of the series 
expansion of the density of states which were scale-invariant. 
Here we demonstrate the scale-invariant regime more gener- 
ally. We define R± by V;ff(i?±) = e (with i?_ = for 
e > Vb - Vm), so that, for Vq > V^, 



2 7->2 



Voe' 



"^'^^ - '^^>l [1 



-ln(yo/K)]. 
(Bl) 



(A7) 



Implicitly differentiating, keeping e and ctq constant, we find 
dR\/dVo = <Tq/Vo is not dependent on e, either directly or 
via R± (which is a function of e). Note that for Vq < Va (the 
dimple or flattened trap), 9i?j./(9Vb is a function of R± so 
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that the following doesn't apply. For e > Vq — Vm, R- — 0, 
so dR^/dVo = in that region. 



From Eq. (37) of 112311 . using 



5(e) 



1 m 



fl4 



^z JR_ 

m 



drr[e- T4ff(r)] 



2^13 



HJ, 



So, 

dg{() ^ _m 
dVo 2hh 



r'(e + KO- j^.V + ^o^o^e-'-^Z-o 



^ (-"-«' -4 



(B2) 

R^ 

(B3) 



^|(4%(^ + ^'"-^o) ifi?-=0 
1 otherwise. 



Hence 5(e) does not vary with Vq for i?_ > 0, i.e. e < Vq 



Vm and the system is scale invariant for /c^T ^ Vq — Vm 
(note that we do not require e ^ Vq ~ Vm)- Interestingly, if 
we define Vq to satisfy Vq = e + Vm, then we can integrate 
( IB4I ) to get the following exact correction to the scale invariant 
density of states, gsi{e) from (l30t , for Vq > Va and e > Vq — 



^n 



^«-^-W 1 (^)3 


(e + K. 


4(>"'|-^f)] 



Vn* 



(B5) 



If we differentiate ^^(F) [Eq. (lA4l i1 



(B4) 



5gr(l") 
5K) 



1 a /(9i?^ 



2 9y V 914 



dRl 

'dVQ 



0, 



(B6) 



then use ( IA2l i. we find that the trap density of states, gtrapiV), 
is also scale invariant for V < Vq — Vm- 
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